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Abstract 

The properties of scattering phases in quantum dots are analyzed with the 
help of lattice models. We first derive the expressions relating the different 
scattering phases and the dot Green functions. We analyze in detail the 
Friedel sum rule and discuss the deviation of the phase of the transmission 
amplitude from the Friedel phase at the zeroes of the transmission. The 
occurrence of such zeroes is related to the parity of the isolated dot levels. 
A statistical analysis of the isolated dot wave-functions reveals the absence 
of significant correlations in the parity for large disorder and the appearance, 
for weak disorder, of certain dot states which are strongly coupled to the 
leads. It is shown that large differences in the coupling to the leads give rise 
to an anomalous charging of the dot levels. A mechanism for the phase lapse 
observed experimentally based on this property is discussed and illustrated 
with model calculations. 



I. INTRODUCTION 

Phase coherence is at the heart of most phenomena studied in mesoscopic physics. How- 
ever, the behavior of the electronic wave-function phase itself in an actual quantum transport 
device had not been studied until recent years. In the case of quantum dots, investigations 
have been predominantly restricted to conductance measurements which carry no infor- 
mation on the transmission phase. It was not until the experiments by Yacoby et al. that 
the interest in phase behavior of quantum dots really started. In these experiments a quan- 
tum dot was embedded in one of the arms of an Aharonov-Bohm (AB) ring in an attempt 
to analyze the transmission phase evolution as a function of the dot gate voltage. Although 
these experiments were the first to demonstrate the presence of a coherent component in 
the current through a quantum dot in the Coulomb blockade regime, they failed to give the 
complete evolution of the phase. This limitation was explained |p as a consequence of the 
phase locking (or phase rigidity) that occurs in a two terminal geometry. It was shown 
that the AB effect in such geometry is characterized by a parity: as a function of the AB 
flux the conductance exhibits either a local maximum at zero flux (positive parity) or a local 
minimum (negative parity). Another intriguing feature of the experimental results was the 



"parity conservation" over a large sequence of Coulomb blockade peaks, which reflected a 
similar evolution of the phase over each peak. The complete evolution of the phase could 
be obtained in a subsequent experiment by Schuster et al. |^ using a four terminal geome- 
try. This experiment confirmed the expected evolution of the phase around the peaks and 
revealed that an abrupt jump of tt occurs in the valleys between the peaks. 

Since 1995 several theoretical efforts have been devoted to explain these observations 
PJ5|-|8|JT(]|,|5|JTT|-p^ : Ref. proposed a screening effect, Ref. alluded to dot degeneracies, 
Refs. associated the observed effect to an asymmetric deformation of the dot which 

leads to repeated charging of the same dot level and Ref. |]TB|] pointed out some special 
properties of the dot states in a semi-chaotic situation (this mechanism will be further 
analyzed in the present work). In spite of all these efforts there is still the feeling that a more 
fundamental explanation is lacking. Each one of the proposed mechanisms can be criticized 
as relying on some particular assumptions. Only the approximate sum rule proposed recently 
in Ref. is supposed to be valid in a generic (chaotic) situation. However, there is still 
no experimental evidence of the near resonance phase-lapse predicted by this mechanism. 

On the other hand, the phase problem affects our knowledge on generic properties of scat- 
tering phases. One fundamental relation, invoked in Ref. in connection to this problem, 
is the Friedel sum rule which relates the phase of the eigenvalues of the scattering matrix to 
the charge accumulated in the dot region. Being related to the dot charge, the Friedel phase 
is a continuous function of the system parameters and cannot exhibit an abrupt behavior 
like the one found in the experiments by Schuster et al. [|]. However, as pointed out recently 
by Lee |jl8| and Taniguchi and one of the authors [0 the phase of the transmission ampli- 
tude can depart from the Friedel phase and exhibit a non-analytic behavior at the points 
where the modulus of the transmission vanishes. It is thus intersting to study the general 
conditions for the occurrence of zeroes in the transmission through a quantum dot. 

The aim of this paper is to investigate the behavior of the different scattering phases 
in quantum dots with the help of lattice models. These type of models allow to describe a 
dot of arbitrary shape and to study the influence of disorder pO[. We shall first derive the 



expressions for the different scattering phases in terms of Green functions. These expressions 
allow to relate the occurrence of zeroes in the transmission with the parity of the isolated 
dot wave-functions. We also study the statistical properties of the dot wave-functions in 
a disordered quantum dot, showing explicitly the absence of significant correlations in the 
chaotic case. On the other hand, as suggested in Ref. , for weak disorder one can identify 
certain dot levels which are much more strongly coupled to the leads than average. We shall 
show that in this situation the dot levels are populated in an anomalous way as a function 
of gate voltage. We shall finally discuss the possible role of this type of correlation effects 
in the phase problem. 

The outline of the paper will be the following: In section II we introduce a generic 
lattice model for a Quantum Dot coupled to single-moded leads. In section HI we derive 
the expressions for the different scattering phases in terms of Green functions. We discuss 
in particular the Friedel sum rule and the relation of the Friedel phase to the phase of the 
transmission amplitude. In section IV we study the conditions for the occurrence of zeroes in 
the transmission and show that they are independent of the strength of the coupling to the 
leads. The statistical properties of the isolated dot wave-functions are analyzed in section 
V. Finally, in section VI we discuss the role of electron correlation effects. We end the paper 
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with some conclusions and final remarks. 



II. GENERIC LATTICE MODEL 

As a model for a two dimensional quantum dot we consider a collection of N sites on a 
square lattice (see Fig. 1). This model can represent a dot of arbitrary shape. The electrons 
in the dot are described by a tight-binding Hamiltonian with site energies and a constant 
hopping element t coupling nearest-neighbors only, {t will be taken as the unit of energy). 
The site energies are allowed to vary following an imposed electrostatic confining potential 
and/or randomly in order to study the influence of disorder. 

On the other hand, electron-electron interactions can be included within the constant 
charging energy model by adding a term Vcoui = Ec^Ndot — CVgjef' , where N^ot is the mean 
number of electrons in the dot, to the one-electron Hamiltonian. Its effect will be discussed 
in section VI. 

We shall consider that the dot is coupled to electron reservoirs by two one-dimensional 
leads as depicted in Fig. 1. We may assume that the coupling to the left and right leads 
is restricted to two sites labeled by 1 and respectively. The first site on each lead are 
connected to these two sites on the dot by hopping elements ti^ and respectively. As will 
be discussed later, this situation can be easily generalized to the case where the first sites 
on the leads are connected to several sites on the dot (this multiple connection is illustrated 
in Fig. 1 by the dashed lines). 



III. SCATTERING PHASES AND GREEN FUNCTIONS 

The electronic properties of lattice models are conveniently given in terms of Green 
functions. We need to introduce the retarded and advanced Green operators given by 



H, 



dot 



(1) 



where Hdot is the one-electron part of the isolated dot Hamiltonian and S*^'" is the (retarded, 
advanced) self-energy operator describing the coupling of the dot to the leads. In a Fermi- 
liquid like description, the self-energy should also contain terms accounting for electron- 



electron interactions ||21]]. We shall postpone its discussion to section VI and concentrate on 
the one-electron properties of the quantum dot until then. When the coupling to the leads 
is localized at sites 1 and we have 



(2) 



where a = L,R, ji = 1, Jr = N and g^"" are the local Green functions at the semi-infinite 
one-dimensional leads. 

In terms of the Green operator one can express the dot total density of states p{uj) as 



—Ira Tr 

71 



(3) 
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We shall now discuss the connection of this quantity to the scattering matrix S{uj) 
defined in terms of Green functions by means of the generalized Fisher-Lee relations p2[ 



S.A^) = 5^,p - 2i^T^TpG]^^^^{u), (4) 

where we have introduced the tunneling rates = t^Im(yf^. The unitarity of S can be 
readily shown (see appendix A). We start by defining the quantity Op as 

Opiuj) = ImlnDet [u - Hdot - . (5) 

The derivative with respect to the energy of is then given by 



— — = ImTr 

OUJ 



(6) 



Thus, in the case when the energy dependence of the self-energy can be neglected one 
obtains the identity 

This case corresponds to a particular type of leads, having a large density of states, which 
can screen any deviation from charge neutrality induced by the presence of the dot. This 
type of leads can be called "non-polarizable leads" . 

On the other hand, it can be show (see appendix B) that Op can be expressed in terms 
of the scattering matrix S as 

ep{uj) = inDet \S{u;)] (8) 

and thus one obtains a relation between the dot total density of states and the derivative of 
the scattering matrix with respect to the energy 

Then, by integrating this expression up to the Fermi energy, we obtain the generalized 
Friedel sum rule 

Ndot = -0f{Ef). (10) 

We see that is an important scattering phase. As it is related to the dot charge it 
should be a continuous function of the energy. It has also a simple relation to the eigenvalues 
of the scattering matrix. Due to the unitarity of 5* its eingenvalues are of the form e^*^^'^ 
and thus Op = + ^2- 

It should be emphasized that relation ([T0|) holds only for the case of non-polarizable 
leads. For a more general case one should include also the charge induced on the leads in 
the Friedel sum rule. The deviation between the dot total density of states and the derivative 
of the Friedel phase with respect to energy has also been pointed out by Gasparian et al. 



4 



23| who analized the connection between densities of states and the scattering matrix for 



continuous models. In particular Eq. (14) in |23] can be written as 



which coincides with our Eq. (^) provided that we make the approximation 

(12) 



Another scattering phase which is relevant for the interference phenomena observed in 
the experiments is the phase of the transmission amplitude 9t = arg Slr. In some particular 
cases (for instance in a one-dimensional problem) Op and 6t coincide. However, as noticed 
recently by some authors [T^|l^, they are in general different. While 6p is a continuous 



function, 9t may not be defined at certain energies where the transmission vanishes. In 
order to be more precise, one can parametrize a general scattering matrix as 

/ ieiie+vi)sm(f) e'^^+'^^^ cos 4> \ , . 

\ e^(''-^2)cos0 ie*(^-^i)sin0 J ' ^ 

with real phases 6, (fi, and 0. It is then easy to show that Op = 9 + ti /2. On the other 
hand, when time reversal symmetry holds one has SpR = Srl and thus ip2 = 0, in which 
case the argument of the transmission amplitude is related to 6*^7 by 

= ^F + 7re(cos0)-^, (14) 

where G(x) is the step function. Therefore, 6t exhibits jumps of vr each time cos0 changes 
sign, i.e. at the points where there is a zero of the transmission. At these points the phase 
of the transmission amplitude deviates from the Friedel phase. Notice that the abrupt 
jump of TT of the phase of the AB oscillations between consecutive resonances is a central 
feature of the experimental results of Schuster et al. 0]. We thus conclude that the study 
of the occurrence of zeroes of the transmission is essential to understand the experimentally 
observed behavior. 



IV. CONDITIONS FOR ZEROES OF THE TRANSMISSION 

Within lattice models one can establish precise conditions for the occurrence of zeroes 
in the transmission amplitude through the dot. According to the Fisher-Lee relations the 
condition for having a zero in SpR at energy Eq is that G[j>^{Eq) = 0. This is gives 

/ 771 \ _ ClwiEo - Hdot - ^''(-E'o)) _ I . 

Det(^o - Hd^t - S'^(^o)) 

where Cij{A) denotes the cofactor of the element i,j in the matrix A. It is easy to see that 
the polynomial in the numerator is real and does not depend on the self-energy coupling of 
the dot to the leads at sites 1 and A^. This is a direct consequence of having the coupling 
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to the leads localized at these sites. We can thus reduce the condition for the zeroes, Eq. 



(15), to the simpler expression 

CiN {Eo - Hdot) = 0. (16) 

Eq. (^) clearly shows that the zeroes of the transmission are characteristic of the isolated 
dot structure and do not depend on the strength of the coupling to the leads. This is 
illustrated in Fig. 2 where the transmission for a 5x5 sites dot is shown for varying values 
of r within an energy range having a zero. One can observe that while the shape of the 
transmission varies substantially, the position of the zero is not affected. 

This property allows one to relate the zeroes of the transmission to the wavefunctions of 
the isolated dot. In the weak coupling limit one can approximate Giat as 

i^'^4-^-A„-zri(^f)2-^r«(^-)2' ^ ^ 

where A„ and denote the eigenvalues and the amplitudes of the corresponding wavefunc- 
tion for the isolated dot. The condition to have a zero between two consecutive eigenvalues 
A„, A„+i is then simply given by > 0. 

We can now identify the sign of ipiip^ as the parity of the corresponding dot wavefunction. 
By this reasoning we conclude that there should be a zero of the transmission in between 
dot states with the same parity regardless of the strength of the coupling to the leads. 

In real systems inelastic scattering would prevent the occurrence of exact zeroes in the 
transmission. This situation can be described within our model by additional leads coupled 
to the dot as voltage probes [^]. This effect is discussed in the next subsection. 



A simple example 

A simple example which already exhibits a zero in the transmission amplitude is the case 
of a four sites sites dot, i.e. = 4. In this model sites 1 and 4 are the ones coupled to the 
leads. Sites 2 and 3 are coupled to sites 1 and 4 by hopping elements t and we take the site 
energies on 1 and 4 as ei = €4 = 0. The transmission amplitude for such model is given by 

Slr = -22^1^^^ ^ ^^^^^^ ^ ^^^^^^ _ _ ^^^^ _ -e){uj + tT) ' ^^^^ 

where e = (es + 63) /2 and f = (E^ + Fk)/2. 

For weak coupling and | £2 — £3 1 ^ t the transmission has two well resolved resonances at 
uj ~ €2,3. The modulus and the phase of Slr are shown in Fig. 3. As can be observed the 
phase of the transmission exhibits a similar evolution around each resonance. At the point 
UJ = e there is a zero of the transmission and its phase exhibits an abrupt jump of vr. From 
Eq. |1^ it is clear that the zero is not dependent on the strength of the coupling to the leads. 

Inelastic scattering can be simulated by additional voltage probes coupled to sites 2 and 
3. Let us denote by Tinei the coupling to the additional leads. It is easy to see that the zero 
of the transmission now moves away from the real energy axis to e + iTinei- The jump in 
the phase of the transmission becomes smaller than vr and it is no longer abrupt but has a 
finite width given by Ti^ei- The width of the phase jump can thus be taken as a measure of 
broadening of the dot levels due to inelastic scattering. 
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V. STATISTICAL ANALYSIS OF THE TRANSMISSION PHASE 



According to the discussion of the previous section, the occurrence of zeroes in the trans- 
mission (and the associated jump in the transmission phase) is related to the parity of the 
isolated dot wavef unctions on consecutive levels. In a similar way as done by many authors 
for analyzing the mesoscopic conductance fluctuations of quantum dots in the Coulomb 



blockade regime |25[, the parity should be studied statistically over many realizations of the 
dot potential. 

In order to search for correlations of the parity of the dot wavefunctions we have nu- 
merically diagonalized dot lattice models up to 30 x 30 sites. The confining potential is 
assumed to have the form of an isotropic parabola with curvature a. We have investigated 
the influence of disorder and also the influence of different models for coupling the dot to 
the leads. 

Figure 4a shows the spacing distribution between consecutive levels for various values of 
the disorder strength W and a = 0.01 (in units of t/a^, where a is the lattice spacing). One 



can notice that the distribution approaches the the Wigner distribution |26 



p[s) = |sexp 



- — s 
4 



(plotted as a full line in Fig. 4a) as the disorder strength increases. This agreement suggests 
that our results for W > 1 should be well described by random matrix theory. 

For analyzing correlations in the parity of the wavefunctions we plot in Fig. 4b the 
probability to find consecutive levels with the same parity as a function of level spacing. As 
can be observed, this probability is almost constant as a function of level spacing. Although 
the probability is slightly larger than 0.5 for weak disorder {W = 0.5) it approaches 0.5 as 
the disorder increases. This behavior indicates that correlations in the parity are negligible 
within this model. 

The previous results correspond to the case where the dot is coupled to the leads at sites 
1 and N, and thus the parity of the wavefunctions is given by aiglipiip^). This situation 
can be generalized to the case where many dot sites are coupled to each lead as indicated by 
the dotted lines in Fig. 1. We shall assume that both leads are coupled to the same number 
Nieads of neighboring sites on each side of the dot. In this case the parity will be given by 
^^siY^jL.ja V'^i^'^^fl)' where Jl^ji denote the sites coupled to the left and right lead respectively. 
It should be emphasized that the conditions for the occurrence of zeroes discussed in section 
IV remains valid in this case, the first sites on the leads L and R playing the role of sites 1 
and A^. 

In Fig. 5 we show the total probability to find consecutive levels with the same parity as 
a function of Nieads and different values of the disorder strength W . For = 0.5 one can 
notice a slight decrease of the probability for increasing Nieads- The probability nevertheless 
remains always close to 0.5 indicating the absence of significant parity correlations even 
when changing the model for coupling the dot to the leads. 

These results demonstrate that it is very unlikely to find a large sequence of dot levels 
with the same parity and consequently the difficulty to account for the experimental results 
of Ref. 1^ within a one-electron model. There are, however, some peculiarities of the dot 
states in the limit of weak disorder which can give rise to correlation effects as discussed in 
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the next section. These correlation effects are associated with certain dot states which are 
strongly coupled to the leads as discussed below. 

The coupling strength = | J^ji^j^ i^j'^i'j'nl shown in figure 6 as a function of the level 
number for the case Ni^ads = 5 and different values of W. In the case of extended contacts 
to the leads and weak disorder one can clearly distinguish states which are more strongly 
coupled to the leads than the average. For the case W = 0.5, this states are indicated by 
arrows in Fig. 6. One can notice by their level number (67,79,92,106,..) that they appear 
in a well defined sequence, corresponding to the shell structure of states in the isotropic 2D 
harmonic oscillator. Contour plots of the wavefunctions with level number 67 and the next 
three levels are shown in Fig. 7. Although they correspond to nearly degenerate levels only 
the first one is strongly coupled to the leads. 



VI. EFFECTS DUE TO COULOMB INTERACTIONS 

Up to now we have neglected interactions and concentrated on the one-electron properties 
of the dot. As commented in section II, electron-electron interactions can be included by 
means of the constant charging energy model. Within the traditional description of Coulomb 
blockade in quantum dots based on simple master equations the main role of the charging 
energy is to open a gap between consecutive dot levels, dot states being populated according 
to the isolated dot ground state for each number of electrons. Within this picture, the 
evolution of the phase over a sequence of dot resonances will be fixed by the parity of the 
corresponding one-electron dot states. As discussed in the previous section, the parity of 
one-electron dot states does not exhibit significant correlations and therefore this picture is 
unable to account for the experimental results of Ref . [|] . 

The situation may change drastically when taking into account the finite coupling to 
the leads beyond lowest order perturbation theory. In this case correlation effects may lead 
to a population of the dot levels different to the one predicted by simple master equations. 



This possibility was recently pointed out by Silvestrov and Imry |16| and advanced as an 
explanation of the experimental results on the phase of the transmission through a quantum 
dot. The mechanism proposed by these authors would take place when there is a strongly 
coupled dot level followed by nearly degenerate weakly coupled levels. This situation is 
characteristic for the case of weak disorder discussed in the previous section. 

The basic mechanism proposed by Silverstrov and Imry can be understood by considering 
a simple four sites problem with two (spinless) electrons. We shall assume that two of the 
sites correspond to the dot levels and the other two represent the left and right leads. Let 
us call ei and e2 the one-electron energies of the isolated dot levels and A = 62 — ei > the 
spacing between them. The lower level is symmetrically coupled to the leads by effective 
hopping elements tn = tm = T which are much larger than the ones corresponding to the 
upper level = t2R = i- The energy levels on the leads sites are taken as zero. Let us 
analyze how the dot levels are populated as the gate voltage shift down the energy levels 
from ei, €2 > to ei, 62 -C 0. 

Finding the ground state of this problem requires in principle to diagonalize the 6x6 
matrix corresponding to the system Hamiltonian in the basis \n1n2nLnfi > where {ni,n2) 
are the occupation numbers for the dot levels and {nL,np() correspond to the leads. The 
problem can, however, be reduced by a basis change in which we replace the states with one 
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electron in the leads by its symmetric and anti-symmetric combination, i.e |1010 >, |1001 >, 
|0110 > and |0101 > are replaced by (|1010 > ±|1001 >)/V2 and (|0110 > ±|0101 >)/V2. 
In this way, the initial 6x6 matrix is reduced into two 3x3 blocks: one block corresponds 
to the empty dot states which couples to the anti-symmetric combinations and the other to 
the doubly occupied dot state coupled to the symmetric combinations, having the form 



/ ei + €2 + \/2T ^/2t 

\f2t £2 

V A/2i ei 



ei v^T 

i/" = I €2 V2t 

a/2T 



In the limit A — > 0, i.e. €\ = —Vg the ground state for each symmetry is given by 



) +2(r2 + t2) 




M +2(T2 + t2) 

When <^ the system starts in the anstisymmetric ground state and the charge of 
the dot levels evolve with Vg according to 



(V^, + A'^)2 + 2(T2 + t2) (y^ + A»)2 + 2(T2 + t2) 

which shows that the dot levels start populating when —Vg ~ T. However, as t ^ T, the 
charge in the weakly coupled level will be negligible as compared to the strongly coupled 
one. 

At = Ec/2 there is a crossing between A" and A"*. As a result the symmetric ground 
state becomes more stable and the charge on each dot level is now given by 

21"^ + iVq + \'Y 2T^ + (Vn + X'f 



{Vg + A^)2 + 2(72 + ^2) ' - (y^ + + 2(T2 + ^2) 

Thus, when crossing Vg = Ec/2 the charge on the strongly coupled state goes from 
~ T'^/{T'^ + t^) to ~ t^/{T'^ + ^2) while the charge on the weakly coupled level does the 
opposite evolution (we neglect small corrections of order (T/EcY). We see that nearly one 
electron is transferred from the strongly coupled level to the weakly coupled one in the 
middle of the charging curve, the total charge in the dot remaining constant (up to order 
{T/EcY). The strongly coupled level is again filled when Vg is close to Ec- The situation is 
similar for finite A <^ T, the abrupt jump is also present but shifted to Vg ~ Ec/2 + A. 

A more realistic description of the actual situation requires the inclusion of many sites to 
represent both the leads and the dot levels. One would need a large number of sites in order 
to simulate a continuous density of states on the leads. Using Lanczos method we are able 
to numerically diagonalize interacting systems with up to 16 sites and 8 electrons. The dot 
levels charge as a function of gate voltage for a systems with 6 sites on a line for representing 
each lead and 4 sites for representing the dot levels, is shown in Fig. 8. The dot levels are 
coupled symmetrically to the leads by hopping parameters ti — 0.1 and ^2 = ^3 = ^4 = 0.02 
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(all energies in units of the charging energy Ec) while the hopping parameter within the 
leads is taken as 0.1. The weakly coupled levels lie at 0.025, 0.030, and 0.035 respectively 
above the strongly coupled level. As can be observed, the strongly coupled level is repeatedly 
charged and discharged in a similar way as for the simple four sites model. One can notice, 
however, that the jumps in the charge of the strongly coupled level are not so much abrupt 
and become progressively less pronounced. 

The next question is how this particular charging of the dot levels affects the transmission 
amplitude and its phase. In order to obtain the transmission through the interacting dot 
we map this finite size cluster into an effective non-interacting system with the same charge 
on each dot level for each value of the gate voltage. In this effective system we use the 
same hopping parameters as in the cluster calculation but replace the leads by infinite one- 
dimensional chains. The effective dot levels are determined self-consistently to get the correct 
level charges. The phase of the transmission amplitude and the transmission probability 
thus obtained are shown in Fig. 9. As expected, the phase follows closely the charge of the 
strongly coupled level. It should be remarked that this behavior is not dependent on the 
parity of the weakly coupled levels. On the other hand, the transmission probability exhibits 
broad peaks that can be associated with the repeated charging of the strongly coupled level 
together with narrow satellite peaks corresponding to the charging of the weakly coupled 
levels. These narrow peaks lie closer to the broader peaks as the cycle number increases. In 
the fourth cycle the narrow peak can hardly be resolved. 



VII. CONCLUSIONS AND FINAL REMARKS 

To summarize, we have analyzed the behavior of scattering phases in quantum dots 
using lattice models. We have first studied the definition of the different scattering phases 
in terms of Green functions. We have shown that abrupt jumps of vr in the phase of the 
transmission amplitude are associated with the occurrence of transmission zeroes. Then, 
we have shown that within lattice models and assuming single moded leads the zeroes are 



independent of the strength of the coupling to the leads p8|. This property allows us to 
relate the occurrence of a zero with the parity of consecutive dot states in the isolated dot. 

We have studied the statistical properties of the isolated dot states as a function of 
the disorder strength. This analysis reveals that there are no significant parity correlations 
between consecutive dot levels. For moderate disorder there appear states which are much 
more strongly coupled to the leads than the average. We have studied charging effects 
in this situation using exact diagonalization of small clusters. We have shown that the 
dot resonances as a function of gate voltage may correspond to the strongly coupled level 
through several cycles. 

The viability of this mechanism as a possible explanation of the "phase problem" de- 
serves further discussion. The presence of strongly coupled levels followed by many weakly 
coupled levels is a reminiscent of the integrable case which gradually disappears for increas- 
ing disorder (see Fig. 6). In this sense, the mechanism is not universal as it would not hold 
for a fully chaotic quantum dot. Thus, this mechanism could be tested experimentally by 
leading the dot into a fully chaotic situation (using, for instance additional gates to distort 
the dot shape, as in Ref. p9[). On the other hand, the transmission as a function of gate 
voltage exhibits narrow satellite peaks close to the broader Coulomb blockade peaks (see 
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Fig. 9) which are not observed experimentally. Being so narrow these peaks could be easily 
washed out by thermal broadening or due to the finite bias voltage used in the measure- 
ments. Of course if certain narrow peaks are not resolved in the experiment, a number of 
additional scenarios might account for parity conservation, even within a non-interacting 
picture. Differing parity conserving scenarios can then still be distinguished according to 
their phase evolution pattern. It is very unlikely that a non-interacting theory can generate 
a a phase evolution like that of Fig. 9 which resembles closely that seen in the experiments. 

One should also comment about the role of spin, not included in the present analysis. At 
very low temperatures spin degeneracy would lead to a Kondo effect. This effect has been 
recently observed experimentally but in somewhat smaller quantum dots and would be 
reflected in the phase behavior as a plateau at ~ 7r/2 in between a pair of resonances corre- 
sponding to the same dot level • At temperatures larger than the Kondo temperature 
spin degeneracy could at most account for the similar phase behavior over two consecutive 
dot resonances. Other spin effects, like Hund's rule require almost perfect symmetry leading 



to orbital degeneracy and so far have been observed only in ultrasmall quantum dots 



To conclude, we point out that our analysis has been restricted to the case of single- 
moded leads. This is an assumption which has been implicitly been used in all the theoretical 
models of the problem following the characterization of the experimental set up given in Ref. 

Nevertheless, the analysis of what would be measured in a multichannel situation is of 
interest and deserves further theoretical investigation. 
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Appendix A 

We show in this appendix the unitarity of the scattering matrix deflned from the Fisher- 
Lee relations. One has 

(S^S)^^ = \1-2iTlG[,\^ + 4TlTr\GI^\'' 

(S^S)^^ = -2zv/r^[(l + 2zTlG1,)GI^ - Gl^il - 2tTnGl,^)] , (19) 
and similar expressions for {S^S)^^ and {S^S)^^ . These elements can be rewritten as 

{S^S) = -22^1^ [GU - GU + 2tTLG'i,GU + 2zr^G?^GW] • (20) 

It is now easy to show that the expressions between brackets vanish identically. This 
fact is a consequence of the identity 
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G1,-G\,-2tTL\G1,\' -2tTn\G 



a I 

in\ 



(21) 



which follows from the definition of the Green operator. Taking the matrix elements (1,1) 
and (1,N) one finds 



a |2 



2iVR\G1p^\ 



Gil ~ — 2iTL\G1i 
Gin ~ G\j^ — 2ir lG1iG\j^ + 2iTjiG1j^G^j^]^ 



(22) 



Appendix B 

We show in this appendix that 2i9F — InDet [5*]. First notice that 



Of = ImlnDet 



u-H. 



dot 



—InDet 

2i 



From the definition of the Green operator one can rewrite this expression as 

Of = ^InDet [1 + (S" - S")G"'] . 



(23) 



(24) 



Finally, using that (E'' — E")jj = —2irL5i^i6j^i — 2irji5i^M5j,N, one can easily check that 



9f = —InDet [S] . 



(25) 
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FIGURES 
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FIG. 1. Schematic representation of a lattice model for a quantum dot. 
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FIG. 2. Transmission probability versus energy for a 5x5 sites dot with Vl = = 0.25 (full 
line), 0.05 (dotted line) and 0.01 (broken line). 
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FIG. 3. Modulus (broken line) and phase (full line) of the transmission amplitude for a four 
sites dot given by Eq. (18) with Tl = Tr = 10 and €2 = —63 = 4. The dotted line shows 
corresponds to the phase when an inelastic tunneling rate Fj^ = 0.1 is introduced. 
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FIG. 4. Level spacing distribution and equal parity probability for a 30 x 30 sites for different 
values of the disorder strength W: 0.5 (triangles), 1.0 (boxes), 2.0 (crosses) and 3.0 (circles). Full 
line corresponds to the Wigner distribution. 
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FIG. 5. Equal parity probability as a function of Ni^ads for W 
(boxes) and W = 2.0 (crosses). 
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FIG. 6. Strength of the coupling to the leads as a function of level number for the case 
^leads = 5 and different vales of W: 0.5 (triangles), 1.0 (boxes) and 2.0 (crosses). The arrows 
indicate states with the strongest coupling to the leads for the W = 0.5 case. 
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FIG. 8. Level charges as function of gate voltage for a 16 sites cluster (see parameters in text). 
The full line corresponds to the strongly coupled level and the dashed, dotted and small dashed 
lines correspond to the consecutive weakly coupled levels. 




12 3 4 



FIG. 9. Phase and modulus of the transmission amplitude for the model of Fig. 8 mapped into 
an effective one-electron problem. 
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